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Abstract. 

As the ground-based gravitational-wave telescopes LIGO, Virgo, and GEO 
600 approach the era of first detections, we review the current knowledge of the 
coalescence rates and the mass and spin distributions of merging neutron-star 
and black- hole binaries. We emphasize the bi-directional connection between 
gravitational-wave astronomy and conventional astrophysics. Astrophysical input 
will make possible informed decisions about optimal detector configurations and 
search techniques. Meanwhile, rate upper limits, detected merger rates, and 
the distribution of masses and spins measured by gravitational-wave searches 
will constrain astrophysical parameters through comparisons with astrophysical 
models. Future developments necessary to the success of gravitational-wave 
astronomy are discussed. 
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1. Introduction 

We are likely only a few years away from the first direct detections of gravitational 
waves from coalescing compact-object (CO) binaries. The Laser Interferometer 



completed two years of data taking at the designed sensitivity (the S5 run) in the 
fall of 2007. Simultaneously, the Virgo detector (http://www.virgo.infn.it, [7]) ended 
its VSR1 run. Other ground-based detectors in operation include GEO 600 [55], 
By 2014-2015, the Advanced versions of the LIGO and Virgo detectors are expected 
to operate at sensitivities ~ 10 - 15 times greater than Initial LIGO [53]. These 
instruments will be able to detect gravitational waves (GWs) from inspirals, mergers, 
and ringdowns associated with coalescences of compact binary systems composed of 
neutron stars (NSs) or black holes (BHs). Although estimates of the rates of such 
coalescences cover a rather wide range, Advanced LIGO may be able to detect tens 
or even hundreds of these events per year (see, e.g., O'Shaughnessy et al. [63]). 

In this paper, we review the mechanisms of formation of coalescing NS and BH 
binaries and summarize the predictions for the merger rates of NS-NS, NS-BH, and 
BH-BH systems. For NS-NS systems these predictions can be made by extrapolating 
from the small sample of observed binary pulsars, while for other binary types they 
are based on population-synthesis models of isolated binary evolution. We discuss 
the potential importance of the dynamical formation of BH-BH binaries in dense 
stellar environments, such as globular and nuclear star clusters, and the possible 
coalescenses involving intermediate- mass black holes (IMBHs). We summarize the 
existing knowledge about the primary characteristics of NSs and BHs in binaries, 
namely, their masses and spins. 

We focus on the tight connection between (i) astrophysical observations and 
modeling and (ii) gravitational-wave searches and data analysis. Input from 
astrophysics can inform decisions about which detector configurations should be 
chosen, which searches are the most likely to yield interesting results, and which 
waveform families should be used for detection. On the other hand, looking beyond 
the excitement of the first detection, the detected merger rate (or even non-trivial 
upper limits) will allow us to limit the space of allowed astrophysical models, thereby 
providing constraints on currently ill-determined astrophysical parameters. Even more 
precise constraints will come from matching the observed parameters of the binaries, 
particularly their masses and spins, with the model predictions. 

2. Mechanisms of Formation and Rate Predictions 

The standard formation channel for compact-object binaries in the field is isolated 
binary evolution. In this case, a primordial binary of main-sequence stars evolves 
through several stages of mass transfer and possibly common-envelope evolution, while 
the binary components age and eventually undergo supernovae, occasionally leaving 
behind a binary composed of neutron stars or black holes (see [U] for a review of this 
process). A circular binary will merge through radiation- reaction from gravitational- 
wave emission on a timescale of [70] 
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where mi and TO2 are the component masses and a is the orbital radius. In this 
section, we describe two different ways to evaluate the merger rates for CO binaries: 
directly from observations for NS-NS systems and via population-synthesis modeling 
for other binary types. We also discuss dynamical formation scenarios that could be 
significant for BH-BH binaries in dense stellar environments. 

2.1. Binary pulsars: observational constraints and extrapolated NS-NS rates 

Among doublc-CO binaries, NS-NS binaries are the only ones for which we have 
direct observational evidence. This evidence comes from observations of Galactic 
binary pulsars: double neutron-star systems where at least one of the NSs is a pulsar. 
Approximately ten such systems have been detected so far. By modeling the selection 
effects of the searches which detected these systems, it is possible to estimate the 
total rate of merging NS-NS binaries in the galaxy, whether or not they are detected 
as binary pulsars (the pioneering efforts to do this are described in [7T] and [ST])- 
Five of the detected systems are expected to merge within a Hubble time through 
radiation reaction from the emission of gravitational waves, and therefore contribute to 
extrapolated merger rate estimates. These systems are: B1913+16 (the Hulse- Taylor 
pulsar), B1534+12, J0737-3039 (a double pulsar), J1756-2251, and J1906+0746; there 
is also a sixth merging system, located in a globular cluster, which is thought to be 
atypical and is not generally included in merger rate calculations. 

Kim et al. [43] developed the statistical framework necessary to reconstruct the 
Galactic NS-NS binary population from observations, while rigorously estimating the 
uncertainties which arise due to the small number of observed systems. Kalogera et 
al. [42] reported the extrapolated rates from the first three systems of the five listed 
above. They used a variety of models for the pulsar luminosity distribution, which 
plays a very significant role in determining the selection effects. The distribution is 
described by two variables: the minimum pulsar luminosity Z/ mm and the slope of the 
luminosity distribution power law p. For the preferred pulsar luminosity distribution 
with -Lmin = 0.3 mJy kpc 2 and p = 2 (see model 6 in Table 1 of [42]), the authors 
found a most likely rate of 83 NS-NS mergers per Myr in the Galaxy; the fifth and 
95th percentile statistical rate estimates were 16.9 and 292.1 per Myr, respectively. 

Although such extrapolation has significantly fewer free parameters than the 
population-synthesis models described below, the numbers should nevertheless be 
used with caution. Different choices of the variables describing the pulsar luminosity 
distribution, which are all consistent with the observations, could shift the merger 
rates by an order of magnitude. The fifth percentile for the model yielding the most 
pessimistic merger rates suggests only 1 merger per Myr in the Galaxy, while rates as 
high as 817.5 per Myr are possible at the 95th percentile for the most optimistic model. 
There are several possibilities for systematic biases related to the imperfect knowledge 
of pulsar age and the beaming fraction }65j . The implicit assumption that the observed 
pulsars form a good representation of the overall Galactic NS-NS population is difficult 
to verify. Lastly, the small number of observed systems means that a single system can 
significantly change the rates predictions. For example, according to [44], the inclusion 
of the relatively short-lived system J1906+0746 could increase the rate estimates above 
by almost a factor of 2. 
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2.2. Population- synthesis models for isolated binaries 

In the absence of direct observations of NS-BH and BH-BH binaries, the only way 
to predict the rates of binary mergers is through population-synthesis modeling. A 
number of different models has been applied over the past two decades (see, e.g., 
[35] 173] |4T| and references therein). Generally, these models rely on a number of 
ill-constrained astrophysical parameters. It is therefore necessary to explore the full 
allowed parameter space in order to quantify the uncertainties in the merger rates 
predictions. We follow the procedure described in [57], since most other models do 
not attempt a thorough exploration of the full parameter space. 

The Star Track population-synthesis code used in |67j has seven free parameters 
that can significantly affect the model outcomes, as identified in [13]. These include: 
the power-law index in the binary mass ratio; 3 parameters used to describe the 
supernova kick velocity distribution; the strength of the massive stellar wind; common- 
envelope efficiency; and the fractional mass loss during non-conservative mass transfer. 
Flat priors are used on these parameters to sample the widest possible set of models. 
However, several additional constraints are then applied to disallow those models (i.e., 
those combinations of the seven input parameters) that disagree with observations. 
The five constraints applied in [67] come from the following observations: (i) merging 
Galactic NS-NS binaries (i.e., binary pulsars expected to merge in a Hubble time); 
(ii) wide NS-NS binaries (those not expected to merge in a Hubble time); (iii) white 
dwarf - NS binaries; (iv) observed Type II supernova rate; and (v) observed Type Ib/c 
supernova rate. Only the first three observations provide non-trivial constraints. The 
predicted Galactic rates for NS-NS coalescences range from 5 to 300 per Myr, with the 
most likely value at 30 per Myr [57]. For NS-BH binaries, the range is from 0.05 to 
100 per Myr, with the most likely value at 3 per Myr [67] . And for BH-BH binaries, 
model predictions range from 0.01 to 30 per Myr, with a most likely value of 0.4 per 
Myr gT]. 

In addition to the ranges stated above, which reflect the impact of different 
choices of astrophysical model parameters, further uncertainties exist. For instance, 
if hypercritical accretion during the common-envelope phase is disallowed, merger 
rates can drop significantly [95] [14] . Many potential BH-BH progenitors may enter 
the common-envelope phase while the donor star is evolving through the Hcrtzsprung 
gap; according to [14] . such systems may merge directly, thus significantly reducing 
the rate of BH-BH mergers. If kicks arc preferentially aligned with pre-supcrnova spin, 
NS-NS coalescence rates could decrease by up to a factor of five [71]. Additionally, 
uncertainties in the star-forming conditions, particularly metallicity, can lead to 
systematic errors in the rates predictions at least at the ~ 30-50% level [68]. Recent 
calculations using updated models of stellar winds and direct observations of the 
black-hole X-ray binary IC10 X-l suggest that low-metallicity environments may be 
even more important [T5][T2]. Meanwhile, all of the applied observational constraints 
suffer from the uncertainties and possible biases discussed in the previous subsection. 
Nonetheless, the default rate predictions for NS-BH and BH-BH merger from all of 
the population-synthesis models mentioned above, as well as those in [58] and [26], 
fall within the uncertainty range listed in the previous paragraph, suggesting that 
predictions generally converge within the uncertainties. 

The rates quoted in Sections 12.11 and 12.21 are computed for our Galaxy. To 
convert them to rates for other galaxies, they can be scaled by the star formation 
rate, which is tracked by the blue-light luminosity for spiral galaxies like the Milky 



Compact Binaries for LIGO and Virgo 



5 



Way [15]. However, while this tracer is relatively accurate for nearby galaxies, it fails 
to account for mergers that arc significantly delayed relative to star formation in older 
elliptical galaxies j^SJ. This omission is particularly significant for BH-BH mergers, 
and therefore the actual BH-BH merger rate in the volume to which ground-based 
detectors will be sensitive in the Advanced LIGO- Virgo era can be significantly higher 
than the contribution from spiral galaxies alone would indicate [63] . 

2.3. Dynamical formation scenarios 

Dynamical effects can make a significant contribution to CO binary mergers in dense 
stellar environments, particularly globular star clusters and nuclear clusters at the 
centers of galaxies. There, two-body encounters can lead to direct captures, while 
three-body single-binary interaction and four-body binary-binary interactions can 
tighten hard binaries to a point where radiation-reaction from gravitational-wave 
emission can take over and drive the binary to coalescence. 

There have been some suggestions regarding the contribution of dynamically- 
formed binaries in globular clusters to short hard gamma ray bursts (SGRBs) [391137] . 
which are believed to be powered by mergers of NS-NS or NS-BH binaries. However, 
extrapolating from the rate of observed SGRBs, as attempted, for example, in [37], is 
highly uncertain, because of unknown selection effects such as SGRB beaming fractions 
[53] . which may be different for field and cluster SGRBs [33]; additionally, binary 
mergers involving NSs are not the only possible formation scenario for SGRBs [94] 146) . 
In fact, most simulations indicate that dynamical effects are not important for binaries 
containing neutron stars, as heavier black holes will segregate to the dense cores and 
substitute in for neutron stars in binaries [79] 00]. In particular, Grindlay et al. [34] 
estimate the NS-NS merger rate from all globular clusters in the Galaxy to be 40 per 
Gyr, a factor of ~ 1000 less than predictions for field mergers. 

On the other hand, dynamical effects can make a very significant contribution 
to BH-BH merger rates, since black holes are the more massive and will sink to the 
central, dense regions of the clusters. O'Leary et al. [61] and Sadowski et al. [79] both 
examined BH-BH binary formation in globular clusters. Their estimates range over 
several orders of ma gnitude, from - 10~ 4 to - 1 BH-BH mergers per globular cluster 
per Myr, due to different assumptions about the fraction of star formation that occurs 
in clusters, the fraction of clusters that may form black-hole subclustcrs, and whether 
these BH subclusters decouple or remain in thermal equilibrium with other stars in the 
cluster. Similarly, Banerjee et al. |11] find that 3-body encounters in intcrmcdiatc- 
age (a few Gyr) massive (> 3 x 10 4 Mq) clusters lead to the formation of a few 
merging BH-BH binaries per cluster, some of which coalesce after ejection from the 
cluster. O'Leary et al. [60] argue that the presence of a massive black hole in galactic 
nuclear clusters can create a sufficiently dense cusp to allow as many as 15 BH-BH 
binaries per Gyr to form in such a cluster through direct 2-body encounters. These 
values may be based on very optimistic assumptions regarding the fraction of BHs in 
galactic nuclei and the number density of low-mass galaxies, but if such BH binaries 
did form, they could be easily distinguished due to their significant eccentricity, which 
is not expected for other stellar-mass binaries in the band of ground-based detectors. 
Meanwhile, Miller and Lauburg [54] calculate that a few x 0.1 BH-BH mergers per 
Myr per galaxy should occur in nuclear clusters of small galaxies without massive 
black holes, primarily through three-body interactions. 

If intermediate-mass black holes (IMBHs) with a mass around 100 Mq are generic 
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in globular clusters (see [33] [55] for reviews on IMBHs, and [55] for an announcement 
of a recently discovered ultra-luminous X-ray source that represents a possible IMBH 
detection), they can contribute to binary merger rates. Since an IMBH will be the 
most massive object in the cluster, it will readily sink to the center and substitute 
into a binary with a compact-object companion. The binary will then harden through 
three-body interactions and eventually merge via an intermediate-mass-ratio inspiral 
on timescales of less than one billion years ,50 . If the binary fraction in a young 
dense cluster is > 10% and the core-collapse time is < 3 Myr, two IMBHs can 
form in a cluster, yielding an IMBH- IMBH merger [55]. Such a merger is also 
possible through the merger of two distinct globular clusters, each containing an 
IMBH [9] . While Advanced LIGO and Virgo will be sensitive to coalescences involving 
IMBHs in globular clusters, third-generation detectors like the Einstein Telescope, 
with improved low- frequency sensitivity |30j . could detect high-redshift mergers of 
IMBHs formed from population III stars which may be the seeds of today's massive 
black holes [E2I331I32]- However, the present lack of knowledge about the existence and 
occupation fraction of IMBHs in globular clusters and about early structure formation 
makes very speculative all rates predictions for mergers involving IMBHs. 

3. Expected characteristics of binaries 

The astrophysical characteristics of merging binaries, namely, their mass and spin 
distributions, are of great interest to gravitational-wave astrophysicists. However, 
these are not well understood. As with merger rates, there are two ways to constrain 
the binary parameters: through direct observations, and through modeling. As 
discussed above, direct observations are limited to a few systems for NS-NS binaries, 
and non-existent for NS-BH and BH-BH binaries. And although BH mass and spin 
measurements exist from mass-transfering X-ray binaries, these do not necessarily 
have the same mass and spin distributions as BHs in compact-object binaries because 
of the particular evolutionary histories of double-CO binaries. Meanwhile, models are 
sensitive to assumptions and can differ significantly in their predictions. We discuss 
what is presently known about mass and spin distributions for NSs and BHs below. 

3.1. NS and BH masses 

The masses of ~ 50 pulsars in binaries have been estimated to date [57] 25J • Until 
a decade ago, most observed binary pulsar masses were extremely consistent with 
a "canonical" value of 1.35M© [55] ■ Since then, several significantly more massive 
pulsars have been found with much larger confirmed masses, such as J1903+0327 
with a confirmed timing mass 1.74M© |24j . In globular clusters, there are four pulsars 
in binaries with weakly constrained orbital inclinations whose masses appear to be 
higher than "canonical", from 1.7M© for a pulsar in Terzan 5 up to 2.74 ± 0.22M© 
for J1748-2021B; see Figure 28 in [47]. 

More than 40 black-hole X-ray binaries have been detected to date. The masses of 
approximately 20 of them have been measured by observing the companion dynamics 
[76] [62]. These masses vary between ~ 4 and ~ 20 solar masses. However, because of 
the mass-transfer episodes involving common envelopes, which play a crucial role in 
the formation of NS-BH and BH-BH binaries, these masses do not necessarily reflect 
the typical BH masses in such binaries. 
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In the absence of direct evidence, masses can be predicted by population-synthesis 
models. Most studies in the literature focus on predictions of the chirp mass, 
M c = (rn 1 TO 2 ) 3 ^ 5 /( m i + "i2) 1 ^ 5 I since this parameter determines the gravitational- 
wave amplitude. The chirp mass is sa 0.87 of the total mass for for equal-mass binaries. 
It is typically found to be between ~ 2 and ~ 5 solar masses for NS-BH systems 
and between ~ 5 and ~ 10 solar masses for BH-BH systems [THJ [66] HH EH [63] . 
However, mean chirp mass values vary by > 30% between these predictions depending 
on assumptions about the cutoff NS mass, accretion rate, etc. The mass ratio is 
typically very poorly constrained in existing population-synthesis studies (see, e.g., 
Fig. 4 of [67]). 

Meanwhile, a number of channels have been suggested through which significantly 
more massive black-hole binaries could form. Recent observations of IC10 X-l, a 
binary with a massive black-hole accreting from a Wolf-Rayet companion star, suggest 
that more massive BH-BH binaries can form through isolated binary evolution, with 
chirp masses ~ 20 Mq [3U] . Dynamical formation channels will favor relatively massive 
BH-BH binary components, as the most massive black holes are likely to sink to 
the centers of clusters through mass segregation and substitute into binaries during 
three-body encounters. Moreover, the BH merger products in such dense clusters 
can be reused if they arc not ejected from the cluster due to recoil kicks, leading 
to higher-mass mergers in subsequent generations; chirp masses for BH-BH mergers 
in globular clusters can therefore range to ~ 30 M Q [79]. Additionally, although 
stellar winds in high-metallicity environments may prevent the formation of massive 
black holes, mass loss through stellar winds would be much less significant in low- 
metallicity environments, allowing more massive black holes to form |12j . Finally, 
IMBHs discussed above could range to hundreds or thousands of solar masses, although 
the more massive IMBHs will not be detectable with ground-based detectors. 

3.2. NS and BH spins 

The observed pulsar spin periods and assumptions about NS spindown rates place 
the NS spin periods at birth in the range 10 — 140 ms [47], which corresponds to 
dimensionlcss spins a/M < 0.04. Meanwhile, the fastest known recycled pulsar in a 
double NS binary has a spin of a/M = 0.02 [23]. Thus, the spins of neutron stars are 
likely to be irrelevant for gravitational-wave detection and parameter estimation. 

On the other hand, direct measurements of BH spins in X-ray binaries through 
quasi-periodic oscillations [86], continuum spectra of the inner edge of the accretion 
disk [SB] , or iron-line profiles [SJ indicate that a significant fraction of BHs are rapidly 
spinning with a/M > 0.7, with some apparently very near to maximum spin, although 
a few could have low spins a/M ~ 0.1. However, all of three methods for spin 
measurement potentially suffer from mo del- dependent assumptions. Moreover, these 
spins could have been enhanced by accretion, and may therefore not be relevant for 
birth spins in NS-BH or BH-BH binaries. 

Population-synthesis models indicate that a majority of both NS-BH and BH-BH 
systems undergo a common-envelope phase with hypercritical accretion, i.e., highly 
super-Eddington accretion with neutrino cooling |41] (however, see |14|). The details 
of hypercritical accretion, including the associated spin-up of the black hole, are poorly 
understood. Under the assumption that all of the angular momentum at the last stable 
orbit is transferred to the BH, which maximizes the BH spin, O'Shaughnessy et al. [M] 
found typical BH spins in access of 0.7, although spins tended to get smaller for more 



Compact Binaries for LIGO and Virgo 



8 



massive BHs, which are more difficult to spin up. However, using more conservative 
assumptions about accretion and spin up, Belczynski et al. |15) conclude that unless 
a BH has a significant birth spin, it is unlikely to be spinning above a/M ~ 0.5. 

The processes governing spin orientations are not fully understood at present. 
For some isolated binaries, the spins may be preferentially aligned with each other 
and with the orbital angular momentum. This will happen if all three of the following 
conditions are satisfied: (i) if the binary was aligned before the supernovae (e.g., 
through mass transfer episodes) , (ii) if the supernovae kicks are small enough that the 
post-kick orbital plane is not too different from the pre-kick one (this is particularly 
likely for the more massive BHs, which are likely to have lower kicks, or for binaries 
that already have very short orbital periods before the second supernova [20], and 
(iii) if the asymmetric supernovae themselves don't change the spin direction relative 
to the spin of the progenitor star. The conditions may be slightly relaxed for the 
primary, as its spin can be re-aligned with the orbital plane through accretion and the 
Bardecn-Pcterson effect. For dynamically formed binaries, there is no clear mechanism 
for preferential spin orientations, so the spin directions should be isotropic. 

4. Gravitational-wave observations and astrophysics 

4-1. Why do we need astrophysics to inform GW searches? 

Astrophysical inputs can significantly aid searches for gravitational waves. Perhaps 
the most obvious such input has to do with the selection of the optimal configurations 
for advanced detectors. Advanced LIGO will have several possibly tunable parameters, 
including the laser power, the reflectivity of the signal recycling mirror and particularly 
the detuning of the signal recycling cavity. Choices of these parameters can 
dramatically change the shape of the noise spectra, as shown in the Advanced LIGO 
Reference Design [3]. For example, the figure in [1] shows six different possible 
configurations. Some of them trade off high-frequency sensitivity for improved 
sensitivity at low frequencies, or broad-band sensitivity for a lower noise floor in a 
narrow frequency band. Future detectors, such as the Einstein Telescope [2], currently 
have an even greater range of possible configurations. Predictions for the rates of 
different classes of compact-binary coalescences and the astrophysical significance 
of specific searches will be very important in selecting detector configurations, 
particularly in the era after the first detections. 

Similarly, astrophysical input can aid decisions regarding the allocation of limited 
human and computational resources in data analysis: Should searches be expanded 
to the mass range of intermediate-mass black holes? Given scarce resources, does it 
make sense to develop a new search for eccentric binaries at the expense of possible 
minor improvements in searches for circular binaries? And what waveforms should 
be used for existing searches [5J [SJ [50]? For instance, the relatively low masses 
predicted for most NS-BH and BH-BH mergers suggest that inspiral-only waveforms 
are sufficient for such searches, but templates that include merger and ringdown 
signals j22j 8] will be necessary to search for mergers involving IMBHs. Although 
current searches do not include spins because existing non-physical spinning template 
families increase the false alarm rates unacceptably [87], the likelihood of significant 
BH spins requires spinning templates to be used for parameter estimation j90) . And 
the computational difficulty of including spinning waveforms in the template banks 
being used for detection may be partially circumvented by ignoring those regions of 
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the parameter space which arc astrophysically unlikely, such as NS-BH binaries with 
very rapidly spinning (x > 0.9) massive (M > 10M Q ) black holes [591IM] . 

4-2. What can GW searches teach us about astrophysics? 

Once the first detections of gravitational waves are achieved, GW astronomy will allow 
us to utilize this new observational window to study the astrophysics of stellar and 
binary evolution and cluster dynamics. 

Some individual detections will already carry very exciting astrophysical 
information. For instance, a single detection with an accurate BH mass measurement 
of > 100 Mq will provide an unambiguous proof of the existence of an intermediate- 
mass black hole. Any measurement indicating a significant eccentricity in a coalescing 
binary would indicate that dynamical effects played an important role in the binary's 
formation. On the other hand, statistics from multiple observations can be compared 
with models to constrain astrophysical parameters. 

Perhaps the most obvious statistic is the observed rate of binary coalescences. As 
discussed in [68], N detections will constrain the rate of mergers to within a factor 
l/yiV, according to Poisson statistics. A Bayesian framework can be used to convert 
observed rates into additional constraints on the population-synthesis models discussed 
in Sec. 12.21 thereby making it possible to more accurately constrain the underlying 
astrophysical parameters. For a given choice of astrophysical model parameters O, 
population synthesis codes coupled to information about galaxy distributions and 
detector sensitivity provide a distribution of the detectable event rate, p(R\Q). If an 
actual merger rate R is measured for a given binary type, then the likelihood that a 
model with parameters O fits the data is 

p(R\0) = J p(R\@) exp (- |i? 2 ~fH dR, (2) 

where ctr includes both the statistical uncertainty in the measurement (equal to 

V^), and the systematic uncertainty in the models. Then the constrained posterior 
probability density function (PDF) of astrophysical parameters is given by the Bayes' 
rule, 

mh) _ £(We), (3) 

p(R) 

where the prior p(O) represents any a priori knowledge about the parameters, and 
p(R) is the evidence which is set by the overall normalization requirement. 

Even in the absence of detections, stringent upper limits can significantly inform 
our understanding of the astrophysical model parameters. Again, a Bayesian approach 
can be used. However, for illustration purposes, we can consider a simple example, 
where the upper limit is employed as a firm cutoff. Consider, for example, the 
population-synthesis prediction for the BH-BH merger rate per Milky Way Equivalent 
Galaxy (MWEG), constrained by the extrapolated observations of merging binary 
pulsars [BB|; this rate is plotted in blue in the left-hand panel of Figure [TJ If 
observations placed a firm upper limit of 1 BH-BH merger per MWEG per Myr, models 
predicting higher rates are ruled out, resulting in the shifted distribution shown in red. 
The right-hand panel of Figured] shows the distribution of the dimensionless parameter 
describing the strength of the stellar wind. Its prior distribution is assumed to be 
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uniform in [0, 1]. The blue histogram indicates that binary-pulsar observations do not 
constrain this astrophysical parameter, as the distribution remains flat (the drops near 
and 1 are features of the plotting routine). However, as the red histogram indicates, 
the additional constraint from the putative upper limit on the BH-BH merger rate 
would preferentially rule out models with low values of the stellar-wind strength, 
thereby significantly informing our knowledge about this astrophysical parameter. 
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Figure 1. The left panel shows the histogram of the BH-BH merger rate per 
MWEG, while the right panel shows the posterior probability distribution function 
of the dimensionless parameter describing the stellar- wind strength (see |66| for 
details). The blue distributions result from the models analyzed in [66] ■ with only 
the observational constraint from coalescing Galactic binary pulsars imposed. The 
red distributions include the additional constraint from a putative hard upper 
limit of 1 BH-BH merger per Myr per MWEG from GW searches. 



In addition to rates, the distributions of the measured parameters of the coalescing 
binaries, particularly masses and spins, can be matched to the model predictions. This 
approach requires several steps. First, for each detected signal, the binary parameters 
have to be measured. This involves measuring not only the best-fit parameters, 
but also estimating the accuracy of the measurement. It is best accomplished by 
sampling the posterior PDF in the parameter space with a Markov Chain Monte Carlo 
[77] 189] [90] or an alternative Bayesian sampling technique [91] ; such techniques have 
been demonstrated to work for the full 15-dimensional parameter space of spinning 
LIGO waveforms [75] . Next, the PDFs from multiple detections have to be combined 
into a single statement about the parameter distribution in the underlying population 
of CO binaries [49], while taking into account selection biases, such as the mass 
dependence of the detection range [16] . At this stage, comparisons can be made 
to the distributions predicted by population-synthesis models. 

The fully Bayesian machinery described above is essential to exploit the 
information available in most parameter measurements, as the low-amplitude 
detections likely in second-generation detectors will not tightly confine their 
permissible ranges. The chirp mass, however, is measured more precisely than 
the scatter expected between independent samples drawn from the (detected) mass 
distribution. In this exceptional case the number of samples, rather than individual 
measurements' accuracy, limits our ability to reconstruct the intrinsic distribution. 
Assuming we do not attempt to exploit other degrees of freedom, the sampled chirp 
mass distribution can therefore be directly compared with the simulated distribution; 
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Figure 2. Using the number of detections only (thin line) and the number 
of detections together with the observed chirp mass distribution (points) to 
distinguish among 238 binary evolution models. The horizontal axis indicates the 
average number fi of binary detections that each model predicts. The vertical axis 
is the number N matc j l of binary evolution simulations (out of the 238 considered) 
whose results are statistically indistinguishable from the reference model, with 
90% confidence. For comparison, the dotted line shows 238/-y//I. 

see for example |21) . A preliminary self-consistent treatment of both the number of 
detections and their chirp mass distribution is shown in Figure [21 

Finally, a wealth of information can be gained from simultaneous gravitational 
and electromagnetic observations of binary coalescences, such as the observations of 
gravitational waves associated with short GRBs [3] . We refer the reader to [7U |B3] and 
references therein for detailed discussions of the scientific promise of multi-messenger 
astronomy, which spans the range from exposing GRB progenitors to enabling precise 
cosmography .59;. 

5. Future developments 

There is a great deal of ongoing work to improve the detection and parameter- 
estimation pipelines within the LIGO- Virgo collaboration. These efforts range from 
the inclusion of new waveform templates and detection statistics in the pipeline |80j 
to improved sampling techniques for parameter estimation [88( I92j to simulated runs 
on numerical-relativity injections as part of the NINJA project [TU] . 

In order to take full advantage of the data that will soon come from gravitational- 
wave detectors, more work is also needed on the astrophysical side. Larger population- 
synthesis runs with improved models are needed to enable comparisons against 
observational results. Most simulations of dynamical-formation scenarios relevant to 
ground-based detectors have not yet quantitatively explored the uncertainties in the 
predicted rates. In general, detailed studies of all input uncertainties in the models 
are needed. 
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We need to prepare for the era of routine detections by developing a framework 
for converting observed LIGO- Virgo merger rates or upper limits for compact-binary 
sources into improved measurements of astrophysical parameters. In order to utilize 
the full wealth of information that will be available in the measured binary parameters, 
such as masses and spins, we must also be ready to combine the parameter PDFs from 
individual events into statements about the population distribution while carefully 
taking selection biases into account, and to compare these observed distributions 
against those predicted by astrophysical models. 

Another direction of theoretical work relates to probing strong-field gravity, and 
perhaps even testing general relativity, with gravitational-wave observations. There 
is a growing body of work on LISA extreme-mass-ratio inspirals as probes of the 
spacetimcs of massive compact bodies (see, e.g., [751 G2J [HSj)- However, it has 
been suggested that even ground-based detectors can be useful tools for studying 
possible deviations from the Kerr metric. For example, gravitational waves from 
intcrmcdiate-mass-ratio inspirals (IMRIs) of NSs and BHs into IMBHs may encode 
enough information to allow comparisons of the mass-quadrupolc moment of the IMBH 
against its expected Kerr value, thereby measuring the "bumpiness" of the black hole 
[17] . Further development effort is needed to make such measurements part of the 
standard data analysis of detected signals. Finally, there are ongoing efforts to study 
the ability of ground-based detectors to probe the equation of state of neutron stars 
through observations of NS-BH coalescences [38] . 
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